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(N : 

^ . Abstract. We investigate regularity properties of molecular one- 

I electron densities p near the nuclei. In particular we derive a re- 

. presentation 

l> ; p{x) = e^^^^p{x) 

^ • with an explicit function JF, only depending on the nuclear charges 

^ \ and the positions of the nuclei, such that /i G C^'^(]R'^), i.e., p has 

Q>^ ■ locally essentially bounded second derivatives. An example con- 

T-H , structed using Hydrogenic eigenfunctions shows that this regularity 

' result is sharp. For atomic eigenfunctions which are either even or 

, odd with respect to inversion in the origin, we prove that /i is even 

■ C^'"(M'^) for all a £ (0, 1). Placing one nucleus at the origin we 

, study p in polar coordinates x = ruj and investigate ■^p{r,uj) and 

I -^p[r,uj) for fixed w as r tends to zero. We prove non-isotropic 

|. cusp conditions of first and second order, which generalize Kato's 

"d ' classical result. 

^ ■ 
> ■ 

■ 

fi: 

> ! 1. Introduction and statement of the results 

^ ' We consider a non-relativistic A^-electron molecule with the nuclei 

' fixed in M . The Hamiltonian describing the system is given by 

N K 
j=l k=l ' ^ ' l<i<j<N ' 

Here the Rk, k = 1, . . . , K, Ri ^ Rj for i ^ j, denote the posi- 
tions of the (fixed) nuclei in with charges Zi,...,Zk, and the 
Xj = {xj^i,Xj^2,Xj^3) G M^, j = 1,...,N, denote the positions of 
the electrons. The Aj, j = 1,...,N, are the associated Laplacians 
so that A = Ylf=i is the 3A^-dimensional Laplacian. Let x = 
(xi, X2, ■ ■ ■ , xn) G IR^^ and V = (Vi, . . . , Vat) denote the points in 



Date: February 7, 2008. 

© 2006 by the authors. This article may be reproduced in its entirety for non- 
commercial purposes. 

1 



2 S. FOURNAIS, M. AND T. HOFFMANN-OSTENHOF, AND T. 0. S0RENSEN 



and the 3A^-dimensional gradient operator respectively. We write 
H = —A + V where V is the multiphcative potential 

^W = -EEr^+ E T^r (1-2) 

j=l k=l ' ^ ' l<i<j<N ' ^' 

Here we neglect the internuclear repulsion W = J2i<k<i<K \r^-Ri\ 
which is just an additive term in the fixed-nuclei approximation. 

The operator H is selfadjoint with operator domain 'D{H) = iy^'^(]R^^) 
and form domain Q{H) = W^^^{R^^) [Hj. 

Let ip be an eigenfunction of H corresponding to an eigenvalue E G 
M, that is, 

Hi; = Ell). (1.3) 

We shall here only consider eigenfunctions ip G L^(M^^)^. 

Note that physical molecular eigenfunctions have to satisfy the Pauli 
principle. This is however irrelevant to our results and we impose no 
such condition. 

The operator in p.lj) (possibly with the addition of the internuclear 
repulsion W) can be considered as the standard model for atoms and 
molecules in quantum mechanics. The analysis of H is fundamental 
for the understanding of the properties of atoms, molecules or, more 
generally, of matter^. 

It is well known that every eigenfunction oi H is Lipschitz-conti- 
nuous PS] and real analytic away from the points in configuration space 
M'^^ where the potential V defined in (ll.2j) is singular (see Section 
7.5, pp. 177-180]). In this paper we investigate regularity properties 
of the electron density p associated to an eigenfunction ip. The density 
p is defined by 

N ^ r 

p{x) = $^P.(a;) = Y. |V'(a:,x,)|'rfx, (1.4) 

where we use the notation 

Xj . (-^l ) ■ ■ ■ ) •^J — 1 ! "^i+l ) • • • ; -^Af ) ; 

and 

dxj := dxi . . . dxj^idxj^i . . . dx^, 



^For scattering problems and for solid state physics, solutions to which are 
not in are also important. 

*^For some problems it is of course necessary to include nuclear motion, and in 
the presence of heavy nuclei relativistic effects have to be accounted for. 
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and, by abuse of notation, we identify (xi, . . . , Xj-i, x, xj+i, . . . , xn) 
with (x, Xj). 

We assume throughout when studying p stemming from some eigen- 
function ip that 

|^(x)| < Coe'^«l^l for all x G R^^ (1.5) 

for some Co, 70 > 0. By [101 Theorem 1.2] (see also [121 Remark 1.7]) 
this implies the existence of constants Ci, 71 > such that 

I VV^(x)| < Ci e-^il'^l for almost all x G M^^. (1.6) 

Since ip is continuous, ()1.5j) is only an assumption on the behaviour at 
infinity. For references on the exponential decay of eigenfunctions, see 
e.g. Agmon Froese and Herbst |Hj, and Simon ^Hl- The proofs of 
our results rely (if not indicated otherwise) on some kind of decay-rate 
for ip] exponential decay is not essential, but assumed for convenience. 
Note that (jl.Sp and p.6|) imply that p is Lipschitz continuous in by 
Lebesgue's theorem on dominated convergence. 

In we showed that p is real analytic away from the nuclei (p G 
C^{R^ \ {Ri,...,Rk}))] for earlier results see also [JU], [H and [Hj. 
Note that the proof of the analyticity does not require any decay of 
ip (apart from G W'^''^{M.^'^)). That p itself is not analytic in all of 
is already clear for the ground state of the Hydrogen atom {N = 
K = l;Ri = 0,Zi = 1): ip(x) = e"'^'/^ so that the associated p 
(up to a normalization constant) equals e~'^'; hence p is just Lipschitz 
continuous near the origin. 

For the atomic case {K = l; Ri = 0, Zi = Z) a quantity studied 
earlier is the spherical average of p which, in polar coordinates x = ruj 
with r = |x| and uj = x/\x\, is defined by 

p{r) = / p{ruj)duj, rG[0, 00). (1.7) 

The above mentioned analyticity result implies that p G C"^((0,oo)). 
The existence of p'(0) and the so-called cusp condition 

p'(o) = -zm (1.8) 

follow from a similar result of Kato for ip itself; see also [Uj and 
[TUl Remark 1.13]. The existence of p"(0) and an implicit formula for 
it was proved in [TUl Theorem 1.11]; see (jl.29j) below for the exact 
statement. 

In the present authors generalized the results of Kato for ip con- 
siderably for the Hamiltonian in (jl.lj) . In the present paper we obtain 
results, partly in the spirit of these findings, for the density p. In par- 
ticular, we prove results on the regularity of the density p at the nuclei 
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and derive identities which the first and second radial derivatives of p 
satisfy. These identities can be interpreted as cusp conditions (analo- 
gously to ()1.8p ). The methods developed in jTU] play an essential role 
in the proofs of these results. 

We indicate the importance of the electron density in quantum me- 
chanics. From the eigenfunction ip it is, in principle, possible to calcu- 
late the energy, various expectation values, etc.; but ip depends on 3A^ 
variables. Physicists and chemists usually aim at understanding atomic 
and molecular properties by means of the electron density which is just 
a function on M? and can be visualized. The density also has an im- 
mediate probabilistic interpretation. 

In computational chemistry density functional methods are of in- 
creasing importance for calculations of ground state energies of large 
molecules. Thereby the energy is approximated by minimizing a 'den- 
sity functional' which depends nonlinearly and nonlocally upon the 
density. The minimizing function is believed to be a good approxima- 
tion to the density itself. The relationship between most of these func- 
tionals and the full iV-electron Schrodinger equation remains unclear 
though. One exception is of course the archetype density functional 
theory, the Thomas-Fermi theory, which is mathematically and phys- 
ically interesting, and very well understood, see [TH] and [12]. For an 
interesting recent review on various mathematical problems related to 
the many models in computational chemistry, see For some work 
on the density p from a numerical point of view, related to regularity 
questions, see 0. 

Questions concerning the one-electron density p, as defined by (jl.4p , 
pose some challenging mathematical problems. Results as given in the 
present paper contribute to a better understanding of the physics of 
atoms and molecules and in addition should have relevance for compu- 
tational quantum chemistry. 

In the following we use the standard definition and notation for 
Holder continuity and Lipschitz continuity, see e.g. [0]. Let / : D 
r2 ^ M, then / G C^'°'{VL) means, for a = 0, that f is k times contin- 
uously differentiable, for a G (0, 1] that the k-th partial derivatives of 
/ are Holder continuous with exponent a. In the case = 0, we often 
write C"(fi) := when a e (0, 1). 

The main result of the present paper is the following. 



Theorem 1.1. Let ip G L^(R^^) be a molecular or atomic N-electron 
eigenfunction, i.e., ip satisfies ()1.3|) . with associated density p. Define 
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K 



:F{x) = -Y,Zk\x-Rk\- (1.9) 

k=l 

Then 

p{x) = e-^(^V(2;) (1-10) 

with 

/iGCi'^(M^). (1.11) 

This representation is optimal in the following sense: There is no func- 
tion J-" : M.'^ M. depending only on Zi, . . . , Z^, Ri, ■ ■ ■ , Rk, but nei- 
ther on N, p, nor E, with the property that e~'^p is in C^(R'^). 

Furthermore, p admits the following representation: 
There exist Ci, . . . , Ck € and z/ : — R such that 

/i(x) = u{x) + 5^ |x - Rk\'{Ck ■ p ), (1.12) 



k=i 

with 



u eC'^^iW) for alia e (0,1). (1.13) 

Remark 1.2. In the case of atoms {K = 1; Ri = 0,Zi = Z), the 
statement of the theorem reads: There exists C G R^ such that 

p(x) = e-^I^V(a;), M = iy(x) + bp (C ■ ^) (1.14) 

with 

u e C2'°(R^) for all a G (0, 1). (1.15) 

To simplify the exposition, we shall give the proof of Theorem 11.11 only 
in the case of atoms. The proof easily generalizes to the case of several 
nuclei. 

Remark 1.3. It will be evident from the proof that the result (ap- 
propriately reformulated) also holds for each pj seperately (see ()1.4|) ). 
The same is true for the results below. 

Proof of the optimality : We study 'Hydrogenic atoms' (A^ = K = 
1] Ri = 0, Zi = Z) and use the notation (contrary to the rest of the 
paper) x = {xi,X2, xs) G R^, r = \x\. In this case, the operator in (jl.ll) 
reduces to Hz = —Ax — Z/\x\. We will present an example where, no 
matter what the choice of (as in the theorem), p = p cannot be 
C^. The argument resembles the proof of the corresponding result in 

El- 
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The Is eigenfunction is "ipisi^) = e~^'''^'^ with Hztpis = ~(^^/4)V'is 
and the associated density is pis{x) = e"^'' . The 2s and 2p eigenfunc- 
tions are 

= (1 - f r)e-^'^/^ ^P^,{x) = x,e-^^'\ 

Both satisfy Hzip = Eip with E = — Z^/16. The associated densities 
are 

P2.(x) = ^lix) = (1 - f r)2e-^''■/^ p,,{x) = ^lix) = x\e-^^l\ 
Consider now ■?/'mixed = + and 

Pmixed = ^mixcd = P2s + P2p + 2V'2^V'2p- 

A simple calculation shows that 

but e^'"pmixed IS just C^'^, siuce the mixed derivative d^J^x^ of 

e^>2s^2p = a;ie^'-/2(l-f) 

does not exist at x = 0. 

But if p = e-^p with p G C^, then 

Pmixed 6 Pmixed Zr 

^ Z 6 Pmixed 

should also be C^, a contradiction. □ 



Note that 'il)2s{x) = 4'2s(y—x) and ip2p{,x) = —ip2p{—x), but their linear 
combination ip is neither even nor odd. 

Remark 1.4. The representation of p as a product p = e'^p with 
a fixed 'universal' JF such that p is by one degree smoother than p 
corresponds to Theorem 1.1 in 1^7^ where a similar result was obtained 
for the eigenfunction if) itself. In that case though, the correponding 
is more complicated since many-particle interactions have to be taken 
into account. For some interesting recent investigation in connection 
with Jastrow factors from a numerical point of view, see [3]. 

The proof of Theorem II. II will be given in the next section. Here we 
just mention that p satisfies an inhomogeneous Schrodinger equation 
whose investigation is crucial for regularity results like the above, as 
well as it was for the results in (TU]. Let H be given by (jl.l|) and 
consider an eigenfunction ip satisfying ()1.3p . To simplify notation we 
assume without loss that ip is real. The equation 

/ V^(x,xj)(iJ-E)^(x,Xj)rfxj = (1.16) 

jR3iV-3 
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leads to an equation (in the sense of distributions) for pj, namely, 

(-^A-f ^)p,+,,^0. (1.17) 

fc=l ' 

Summing (jl.l7p over j we obtain the equation for p, 

k=l ' 

with h = Ylf=i^j- '^^^ functions hj will be given explicitely in Sec- 
tion see (1231) • 

In |10j we considered the spherically averaged density p (as defined 
by ()1.7p ) for the atomic case. The regularity of h (the spherical average 
of h above) was crucial for the results obtained there. Here we study 
the non-averaged density p for the general case of molecules. Again, 
the regularity of h is essential for our results. 

We continue to consider p in the neighbourhood of one nucleus with 
charge Z. Without loss we can place this nucleus at the origin. 

The equations ()1.12j) and ()1.13|) show that it is natural to consider 
the behaviour of p{ruj) for fixed u as r tends to zero. 

Theorem 1.5. Let ip G L^(]R^^) be a molecular or atomic eigenfunc- 
tion, i.e., ip satisfies (jl.3|) . with associated density p. Assume without 
loss that Ri = and write Z instead of Z\. Let tq = minfc>i 
(Vo = oo for atoms) and let G be fixed. 

(i) The function r i— p{r,uj) := piroj), r G [0,ro), satisfies 

p(-, u) G C2'"([0, ro)) for all a G (0, 1). (1.19) 

(ii) Denote by ' the derivative and define 

r,{x) = e^l-lp(x) , X = ^ - r\C ■ u), (1.20) 

where C & is the constant Ci in ()1.12|) (resp. C in ()1.14j) ). 

Then 

7] G Ci'i(5(0, ro)) , X e C2'"(5(0, ro)) for all a G (0, 1), (1.21) 
and 

p'(0,^) = -Zp(0) + uo ■ (Vr7)(0), (1.22) 

p"(0,c^) = Z2p(0) + 2uj-[C- ^(Vr;)(0)] + u ■ {{D\){{})uo) . (1.23) 

Here (D^x) (0) is the Hessian matrix of x evaluated at the ori- 
gin. 

Remark 1.6. 
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(i) For atoms, rj equals /z from Theorem 11.11 and x equals v from 
Remark 11.21 

(ii) Note that p.l9|) trivially implies that p(r, u) = e~^^rj{r, uj) with 
T]{-,uo ) G C2'"([0,ro)) for all a G (0,1). Compare with (fCTlll . 

(cm). 

(iii) In ^01 Theorem 1.11] it was proved that p defined by ()1.7p 
belongs to C2([0, ro))nC2'°((0, tq)) for all a G (0, 1). (The proof 
in ^U] for the atomic case easily generalizes to the molecular 
case.) Reading the proof of ^lOj Theorem 1.11] carefully, one 
sees that it in fact yields p G C2'"([0, Tq)). The statement in 
(jl.l9|) shows that for fixed c<j G §^ this holds already for p(-, tu), 
i.e., without averaging. 

(iv) The identities p.22|) and ()1.23|) can be considered as non- 
isotropic cusp conditions of first and second order. They 
generalize the cusp condition ()1.8|) . as well as the previously 
mentioned result in |10J for p"(0); more on this in Remark 11.81 
(ii) below. See also the second order cusp conditions obtained 
in [7j for the eigenfunction ip itself. 

(v) It is worth noting that (jl.22|) and (jl.23p can be interpreted as a 
structural result for the density p: From Theorem II. 51 it follows 
that in a neighbourhood of a nucleus (which is at the origin), p 
satisfies (for all a G (0, 1)) 

p(r, uj) = p(0) + r0i(cj) + rV2(c^) + 0(r2+") , r i 0, (1.24) 

and ()1.22j) . ()1.23p show that (pi is a linear and 02 a quadratic 
polynomial restricted to S^. 

It is a natural question whether (jl.24jl extends to higher orders. 

We continue with the atomic case. In view of Remark 11.21 p.l4|) 
and the considerations after the proof of the optimal regularity of p in 
Theorem ll.il the following theorem is natural. 

Theorem 1.7. Let ip G ^^(R^^) be an atomic eigenfunction with as- 
sociated density p. Suppose that 

|^(x)| = |V'(-x)| for all X G M^^. (1.25) 

Then p satisfies 

p{x) = e-^l^lp(x), p G C^'^iR^) for all a G (0, 1). (1.26) 
Furthermore, 



p'{0,u) = -Zp{0) , p"(0,cu) = ZV(0) ((DV)(0)cu). (1.27) 
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We also have 

p"(0, u) = I (ZV(0) + MO, ^)) + \ lim i^Mp^^ (1.28) 

with h from p.lSj) . and CT'jr^ the angular part of —A, i.e., A = 
dydr^ + {2/r)d/dr-Cyr\ 

Remark 1.8. 

(i) In this case ^ = v = x = V-i^^ can be seen from Remark 11.61 
(i) and the proof of the theorem. 

(ii) Note that p.27|) shows that the cusp condition p.Sp in this case 
holds for fixed angle G §^ without averaging. Further, taking 
the spherical average of p.28|) . we get the formula for p"(0) 
obtained in |10^ Theorem 1.11 (iv)]: 

P"(0) = ^(^V(0) + M0)). (1.29) 
To see this note that for all r > 

/ I - {C'^p){r,oj)doj = I (Ch) ■ p{r,uj)duj = 0. 

Note that p"(0) > 0, since 

h{r) > ep{r) 

for some e > JOl Theorem 1.11]. This positivity is not an 
obvious consequence of the formula in ()1.27|) . 

(iii) As can be seen from the proof of Theorem 11.71 h G C"(]R^) for 
all a G (0, 1) in this case. 



2. Proofs 

Proof of Theorem M . 1\ : As noted in Remark 11.21 we shall give the proof 
only in the case of atoms {K = 1; i?i = 0, Zi = Z). 

For the regularity questions concerning p defined in fll.4|) it suffices 
to consider the (non-symmetrized) density pi defined by 

Pl(x) = / \lp{x,X2,...,XM)\'^dX2---dxN 

jR3iV-3 

= / |V'(x,xi)|2c/xi (2.1) 

7R3iV-3 

with X G M^ xi = (X2, . . . , x^) G M3JV-3_ 
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As explained in ()l.lfij) - ()1.18j) pi satisfies the Schrodinger-type equa- 
tion 

1Z 

- Api - — pi + 2/ii = 0, (2.2) 

where the function hi is given by 
hi{x) = Ji- h-Epi{x), (2.3) 

N ^ r z 

M^) = y2 \VM^d±i , J2{x) = V / --^^dki, 

~IJr3N-3 ~^JR3N-3\Xj\ 

^ f 1 r 1 

J3(x) = V/ :1p^d±i+ V / r1p^d±i. 

^ J^3N-3 \X-Xk\ 2<f^<N I ~ I 

(We will henceforth partly omit the variables in the integrands). Using 
the exponential decay of (jl.Sp and of Vip ()1.6|) one can prove that 
hi e L°°(M3) (for details, see [ffl Theorem 1.11]). 
Making the Ansatz 

Pi(x) = e-^N^i(x) (2.4) 
and using ()2.2p . we get that pi satisfies the equation 

Api = 2e^l^'l/ii + 2Zuj ■ Vpi - Z^i- (2.5) 

Here uj = Since pi G C"^'^(M'^) (as mentioned in the introduction), 

also Hi G C°'^(M^). Clearly, the function x ^— belongs to L°°(]R^). 
The fact that also x ^— hi{x) is in L°°(]R^) gives, by standard elliptic 
regularity f7| Theorem 10.2], that 

/ii G C^'"(M^) for all a G (0, 1). (2.6) 

Our aim is to prove more, namely that 

Api = ci ■ + ^ , ci G M^ ge C°(M3) for all a G (0, 1). (2.7) 

Since, by (j2.6j) . Vpi is continuous at the origin, the term 2Zuj ■ Vpi 
behaves like c^i^ ■ uj {c^^^ = 2ZV/ii(0) G M^) at the origin. It turns 
out that generally hi is discontinuous at the origin, also behaving like 
cf'^ ■ UJ (cp"* G M^). However, one can solve the equation Aui = ci ■ uj 
(ci G M^) explicitely, and one gets that the solution ui is C^'^. From 
standard elliptic regularity, the other terms give contributions which 
belong to C^'"(]R^). Below we give the details. 
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First, we write 

2Zu ■ Vfii{x) = ■ iu + gi{x) , gi E C°(R^) , « G (0, 1), (2.8) 

cS'^ = 2ZV^i(0) , (7i(x) = 2Zuj ■ (V/ii(x) - V^i(O)). 

That gi G C"(M^), a G (0, 1), follows from Lemma 1X11 in AppendixEl 
We next consider hi as defined in ()2.3|) . We will show the following: 

Lemma 2.1. Let hi be as in Q. Then there exist cf^ G M^ G : 
]R3 R, such that 

hi = cf^ -uj + G, Ge G°(M^) , a G (0, 1). (2.9) 

Before proving Lemma 12.11 we finish the proof of Theorem 11.11 
Lemma f2. II and Lemma fA. II in Appendix 1X1 imply that 

(e^l^'l - l)hi G C°(M^) for all a G (0, 1), 

and it therefore follows from (ESI), (EH), and (0 that 

Afii = ci ■ u + g , ci = cS^^ + 2cf ^ G , (2.10) 

^ G C"(M3), a G (0,1). 

This is (12. 7|) . which we aimed to prove. 

A simple computation shows that the function ui{x) = ||xpci ■ u = 

^ ci ■ X satisfies Ami = ci ■ tu, and so ui = fii — ui solves Az/i = g, 
g G C"(M^), a G (0,1). From standard elliptic regularity theory [T7| 
Theorem 10.3] follows that ui G C2'°(m3), a G (0, 1). Note that due to 
Lemma Em ui G C^'^(]R^), and so 

/ii = z/i + Ml = i/i + l\x\\ci ■ u) G C^'^(M=^). (2.11) 
o 

This finishes the proof of Theorem 11.11 for atoms, with 

1 ^ 

C = -5^c„ (2.12) 

i=i 

where Cj is the contribution from pj. □ 
It remains to prove Lemma 12.11 

Proof of Lemma \2.1\ : The proof is essentially a tedious but elementary 
verification, the idea being to isolate and extract the most singular 
term of hi. Part of this has been carried out in ^U], and, in order not 
to repeat the details, we refer to that paper whenever possible. We 
also use the same notation. 
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Define 

^1 = e-(^-^iV, (2.13) 

with 



^W = -f El^^-l + i E l^.'-^^l' (2-14) 

j=l l<j<k<N 
„ N ^ 

^i(x) = -^5^V^k,f + 1 + - 5^ + (2.15) 



2 ^ V ' ^' 4 

jr = l l<i<fe<iV 

Tlien it follows from [IQ, Proposition 1.5] that ipi G ^^'"(M^), a G (0, 1). 

In proi, Lemma 3.5 (i)] it is proven that (with the notation from ()2.3|) ) 
J25 J3 £ C°(]R'^), a G (0, 1). Furthermore, using ()2.13p . Ji is written as 

J^{x) = I I rfxi = Ji + /2 + /s + /4 + /5 + h, 

Jr3N-3 

where 

AT 



N 

J2 / iV.FilVc^xi, 

~7 Jr3N-3 



1 

N 

Io(x 



3-- 

N 



TV 

N 

2E / (V,F-V,^i)e2(^-^^)^icixi, 

-2 V / (V,Fi ■ V,^i)e2(^-^^)^i rfxi. 

~^ JIR3JV-3 



J 

N 

h{x) = 

3 

N 

3 

N 

hix 



It is proven in ^01 p. 93, bottom] that 

/2,/4,/6eC"(R3) , aG(0,l). 
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It is also proven in that 

h{x) = f |ViF|Vrfxi + /i(x), 



h{x) = -2 / (ViF ■ ViFi)^=^ d±i + h{x), 

h{x) = 2 / (ViF ■ ViV^i)e2(^-^i)^i dxi + i,{x), (2.16) 

JIR3JV-3 

with Ij e C"(R3),a G (0,1). (For h, this is [IHl (3.52) and until 
(3.53)]; for I3 and J5, this is [TU[ (3.55) and (3.56), and between]). 
Now, (see (HJ^ ) 

I , , ^lO/ ^ Zj X V ^ X Xj 1 / V ^ X X'i 

|V,FP(x.xO = ---^-Era+16(5: 



4 4 a; x — a;J 16V'^ a; — x 

i=2 ' i=2 I 



Therefore, we write 



/ |ViF|Vc^xi = /i,i(x) + /i,2(x), 

/i,i(^) = -7n- / E^^^'^^i' 

4 X Jk3JV-3 ^ X - Xj 



For III, we write /i,i(x) = ^fi ' ^^'^lix), with 



Jk3N-3 l 4 ISV-^IX — Xj|/ J 

The proof that Ii^2 belongs to C"(R'^), a G (0, 1), is also in [lOj (between 
(3.52) and (3.55)'). 

4 \x\ 

Inti(x) = -/ y ^\ ^^d±i. (2.17) 

The function Inti belongs to C"(M'^), a G (0, 1). This follows by argu- 
ments as for the integral 

1 



■ (x ■ (x — Xfc)-?/'^) dx2 ■ ■ ■ dXN 

|X — Xfcl 

in jHH (see between (3.54) and (3.55)). 
Therefore (with 00 = -^^^), 

/i,i(x) = ci,i • CO + (?i,i(x) , (71,1 G CiR') , a G (0, 1), (2.18) 
ci,i = |lnti(0) , gi,i{x) = ■ (inti(x) - Inti(O)). (2.19) 
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That i G C"(M^), a G (0,1), follows from Lemma lA.ll in Appen- 
dix |X1 

Next we consider J3 (see fl2.16|) ). Since 



7 1 ^ _ 



2 a; 4 ^-^ \x — x 

7 1 ^ _ . 

ViFi(x,xi) = --^=== + - ^ 



we have 



~2 I I ^M3^-3 V F ^ ^il + 1 
42(x)= / {^^+(iy^^).ViFiy^xi. 



i=2 



That the first term in 1^,2 belongs to C"(M^), a G (0,1) follows by 
arguments as in llOj (by applying Lemma 3.4 as done after (3.50); note 
that the function X ^ kl/Vkp + 1 belongs to C0'1(R3)). 

That the last term in 1^,2 belongs to C"(]R^), a G (0, 1), is proved in 
Uni (3.55), and after]. 

For /s,!, we write /3,i(x) = -f |f| ■ Int3(x), with 



N _ XX- 



Int3(x) = ^ / , ^ ^'d±,. (2.21) 

Similar arguments as for the first integral in I3 2 above show that Int3 
belongs to C"(R3), a G (0, 1). 

This implies, by Lemma fA.H that 

h,iix) = C3,i ■ uj + (73,i(x) , (?3,i G C°(M3) , a g (0, 1), (2.22) 
C3,i = -|^Int3(0) , (73,i(x) = -^uj ■ (Intsix) - Int3(0)). (2.23) 
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Finally, we consider (see ()2.1fi|l ). We use the same kind of analysis. 
Using (|2.2(J|) . we write 

/ (ViF • Vi^i)e2(^-^^Vi rfxi = /5,i(x) + 15,2(0;), 

jR3iV-3 

/5,i(^) = -f A - / e^^^-^^)Vi(V^D^xi, 

^ N I, _ 

aY. / ^^■(Vi^i)e^(^-^^Vicix,. 



AT 

1 ' 



4 

Again, the proof that 15,2 belongs to C"(M^), a G (0,1), is in [TUl 
(3.56), and after]. As before, we write /5,i(a;) = ~f |f| " Int5(x) with 

Int5(x)=/ e2(^-^^)Vi(^?)dxi. (2.24) 

JIR3JV-3 

The integral Ints belongs to C"(M'^), a G (0,1). This follows by ar- 
guments as in PU] (the term in (3.56) with j = 1). Therefore, by 
Lemma lA.ll 

h,i{x) = C5,i ■ uj + (75,i(x) , ^?5,i e C"(M3) , a g (0, 1), (2.25) 

C5,i = -|lnt5(0) , g5,i{x) = -^u ■ {inhix) - InUiO)). (2.26) 

It follows from all of the above that 
h = (ci,i - 2c3,i + 2c5,i) ■ 00 + [{gi,i - 2^3,1 + 2^5,1) (2.27) 
+ (/i,2 - 2/3,2 + 2/5,2) + (/2 + h + h) 

+ (/i + /3 + h) -J2 + J3] - Epi 
= (ci,i - 2c3,i + 2c5,i) -uj + f- Ep^, with / G C"(M3) , « g (0, 1). 

Since pi G ^^.^(mS), (jTTTjl shows that /ii indeed can be written as in 
(j2.9|) : that is, this finishes the proof of Lemma f2. 11 □ 



Proof of Theorem li.^ : That rj G C^'^(/?(0, rg)) follows from p.9p - 
firm) . That X e C2,«(5(o,ro)),a G (0, 1), is a consequence of (ITT^ - 
firni) . It remains to prove dTTni) . (irail . (ir^ . 

As for Theorem ll.il we shall only give the proof for the case of atoms 
{K = 1; Ri = 0, Zi = Z,ri = fi and x = z/). 

Let in the sequel = G be arbitrary, but fixed. Note that 
()1.14|) and p.l5|) imply that p(r, u) = e~^^/i(r, u) with p( ■ , cu) G 
C2'°([0,oo)). It follows that p{- ,u) G ^^'"([O, cx))) since r e"^'' 
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belongs to C°°([0, oo)). In particular, p'(0,co') = \imrio p'{r,uj) and 
p"(0, uj) = lim^jo P"{r, uj). (All the above for all a e (0, 1)). 
Next, by the above, 

lim p'(r, u) = lim \u! ■ Vp{r, u)] 

rlO rj,0 

= lim [ - Zp{r, uj) + e'^^'tu ■ V/i(x)] = -Zp{0) + uj ■ Vp(0), 
which is (fr^ . 

Finally, the proof of (fT^ . Due to p = e-^> and (im) - (irTH|) we 
have 

p"(r,cj) = Z'^p{r,u) - ^Ze-^'u ■ Vp(s) 

+ e-^'' {uj ■ V(cj ■ Vz/(x)) + 2^ . cj] . 

A simple computation shows that uj ■ Viyj ■ Vz/(x)) = uj • ((i5^z/)(x)co'), 
and, since p G ^^'^(M^) and i/ G C2'"(m3)^ ^ ^ j^g, 1), we get (IT^Hll . 

□ 



Proof of Theorem |i. 7| : We will show that the symmetry assumption 
flT^ for implies that ci = G MMn ^l^i. Then 

Api=^7, (7GC"(M3),ae(0,l), (2.28) 

and so standard elliptic regularity implies that pi G C^'"(M'^) for all 

a G (0, 1). This will prove (IT^. 

Recall that (see (EIHl), (EH), (EUHl), (EH, and ((21211)) 

Cl = cl') + 2cf ) = cS') + 2ci,i - 4c3,i + 4c5,i, (2.29) 
cS') = 2ZVpi(0) , ci,i = jlnti(O), (2.30) 

C3,i = -|lnt3(0) , C5,i = -f Int5(0). (2.31) 
o 4 

We first consider c'l^ . The assumption p.25|l clearly implies that 
pi is an even function on M^. It follows from pi = e^'^'pi that pi G 
C^'-'^(R^) is even and therefore Vpi G C°'-'^(M^; M^) is odd. In particular, 
Vpi(O) = 0, and so S^^ = 0. 

It was shown in the proof of Theorem II . II that Intj G C"(R^; M^), J = 
1,3,5, a G (0,1) (see (ETTj) . ^Tm . and (j2211))- Furthermore, the 
symmetry condition (jl.25j) clearly implies that all three functions are 
odd (for Ints, use ()2.13p ). It follows that Intj(O) = 0, j = 1,3,5, and 
therefore (see (jTipi and ((2211)) c^-i = 0, j = 1,3,5. Therefore cf^ = 
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and hence ci = in ()2.29j) . This, and ()2.7|1 . imphes ()2.28j] . which, as 
mentioned above, proves (jl.26j) . 

Note that cf^ = imphes that h e C"(R3) for ah a E (0, 1) (see 
Lemma 12.11) . 

The above clearly implies that V/i(0) = and C = Ylf=i — O5 ^^"^ 
so (fT^ and (fT^ i mply (ITTfll . 

It remains to prove ()1.28|) . With A = d"^ /dr"^ + {2/r)d/dr — C? jr'^, 
()2.2|) becomes (after multiplication by — r) 

rp"(r, uj) + 2p'(r, cu) + 2Zp(r, uj) - 2r/i(r, cu) = (^'/^)(^'^) , 

r 

This implies, using the fact that h e L°°(U.^) and (f07|) . that 
^.^ (£2p)(r^ ^ 2(p'(0,cu) + Zp(0,a;)) = 0. 

Let 7^(r,cJ) := then reads 

2 

p"(r, + - {p\r, Lu) + Zp{r, u)) - 2h{r, cu) = n{r, cu) , r > 0. (2.32) 
Note that, by I'Hopital's rule and (fTTTIl . 

lim - (p'(r, uj) + Zp(r, cu)) = 2 (p"(0, uj) + Zp'(0, cj)) , 

riO r 

and so ()2.32p implies that lZ{Q,uj) := lim^ jo "^(^^5 1^) exists, and 

7^(0, cj) = 3p"(0, to) + 2Zp'(0, cu) - 2h{0, to). 

The existence of h{0,to) := lim-^jo ^('", follows from Lemma |2.1l 
Therefore, using (jl.27|) . we obtain (jl.28|) . □ 



Appendix A. A useful lemma 

The following lemma is Lemma 2.9 in [7j we include it, without 
proof, for the convenience of the reader. (The proof is simple, and can 
be found in [7j). 

Lemma A.l. Let G : U ^ for U C M""*""* a neighbourhood of 
a point (0,?/o) G M" x M™. Assume G{0,y) = for all y such that 
{0,y) e U. Let 

j §-,-Gix,y) x^O, 
^^"'^)= " . = 0. 
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Then, for a G (0, 1], 

G e C°-"{U; W) ^ f e C^'^^iU). (A.l) 
Furthermore, ||/||c"{c/) < 2||G'||c•«(c/)■ 
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